Valence-bond-solid order in antiferromagnets with spin-lattice coupling 
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We propose that a valence- bond- solid (VBS) order can be stabilized in certain two-dimensional 
antiferromagnets due to spin-lattice coupling. In contrast to the VBS state of the Affieck-Kennedy- 
Lieb-Tesaki (AKLT) type, the spin-lattice coupling-induced VBS state can occur when 2S is not an 
integer multiple of the coordination number z. S = 2 spins on the triangular lattice with z = 6 is 
discussed as an example. Within the Schwinger boson mean-field theory the ground state is derived 
as the direct product of states, one of which represents the conventional long-range ordered spins, and 
the other given by the v3 x y3 modulation of the valence bond amplitudes. Excitation spectrum 
for the modulated valence bond state is worked out within the single-mode approximation. The 
spectrum offers a new collective mode, distinct from the spin wave excitations of the magnetically 
ordered ground state, and observable in neutron scattering. 

PACS numbers: 75.10.Hk, 75.10.Jm 



Spin-Peierls phenomenon refers to the dimerization of 
antiferromagnetic 5 = 1/2 spins on a linear chain ac- 
comapanied by spontaneous lattice distortion [1, 2]. The 
physics underlying the spin-Peierls transition is the extra 
gain in energy through the reduction of dimerized bond 
length. A similar phenomenon is believed to occur in 
two or higher dimensions [3]. Another route to obtain- 
ing dimerized ground state is by including interactions 
beyond the nearest-neighbor Heisenberg exchange, as in 
the Majumdar- Ghosh model for one-dimensional antifer- 
romagnetic chain[4]. 

Generalization of Majumdar- Ghosh idea to two di- 
mensions and to arbitrary spin 5, and construction of 
a general class of exact ground states to antiferromag- 
netic spin models were carried out by Affleck, Kennedy, 
Lieb, and Tesaki (AKLT) [5]. AKLT's idea, and the sub- 
sequent refinement by Arovas, Auerbach, and Haldane 
(AAH)[6], emphasize the close connection between the 
spin quantum number 25 and the lattice coordination, 
z. A valence-bond-solid (VBS) state is formed when 25 
is an integer multiple of z, and every nearest-neighbor 
bond is covered by 25/ z dimers. 

In this paper, we address the case where the two 
numbers, 25 and z, are incommensurate. An example 
would be 5 = 2 antiferromagnet on a triangular lat- 
tice, where the AKLT condition would require 5 = 3. 
Such a system is realized, for example, in the compounds 
RMn03(R=Lu,Y)[7-9]. In these systems, each magnetic 
Mn ion carries spin 5 = 2. At most four valence bonds 
can be formed from each site, leaving the other two bonds 
"unfulfilled" . Long-range magnetic order has been found 
in these compounds. Furthermore, there is a strong in- 
dication of the movement of the Mn position concurrent 
with the magnetic transition, pointing to the existence of 
spin- lattice coupling in these materials [9]. 

Attempts to understand the effect of spin-lattice cou- 



pling for the triangular network of Heisenberg spin for 
large 5 has been given in Ref. 10. It was shown that the 
mean-field theory predicts lack of lattice displacement as 
the magnetic order builds up, in contrast to the exper- 
imental observation. Other types of order, such as the 
valence bond order, had not been examined in detail. It 
is shown here that a type of VBS order with only four 
of the six bonds of the triangular lattice being filled by 
singlets can be stabilized through spin-lattice interaction 
and leads to lattice deformation as in the spin-Peierls 
systems. 

Our ideas are illustrated in Fig. 1 for 5 = 2 and z = 6 
(triangular lattice). Such states are energetically unfa- 
vorable in comparison to the magnetically ordered state 
which makes use of the exchange energy for all six neigh- 
bors. On the other hand, the spin-lattice coupling can 
reduce the distance for a pair of sites forming a valence 
bond, and increase it for sites lacking a bond. The gain in 
the exchange energy achieved in this way will overcome 
the loss from the two missing bonds for a sufficiently large 




FIG. 1: Realization of partial VBS order on a triangular lat- 
tice with spin S = 2. At most four valence bonds can emanate 
from each site (thick lines), forming a regular array of "miss- 
ing bonds" (dotted triangles). Magneto-elastic coupling can 
alter the bond distance depending on the presence/absence of 
the valence bond. Three inequivalent sites are denoted A, B, 
and C. 
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coupling strength, and the VBS state is energetically fa- 
vored over the magnetically ordered state. 

We provide a quantitative justification of our claim by 
analyzing the lattice-coupled spin Hamiltonian, 

tf = X> s *- s i + f E u ?' (i) 

(ij) 1 

taking = J (d° /djj) 7 with the equilibrium and actual 
ionic separations given by d\j and dij, respectively. For 
transition metal ions 7 is a rather large number, ranging 
between 6- 14 [11]. Displacement of the i-th site is given 
by Ui, which costs an elastic energy of (K/2)uf. For ease 
of analysis, we Taylor-expand Jij in small displacements 
as Jij = Jq[1 — jeji ■ (uj — Ui)], taking d®j = 1, and 
£ji = ( r j ~ r i)/\ r j ~ r i\ 1S the urut vector connecting the 
equilibrium ionic sites, rj and rj. 

In the following we re-define Ui — ► u; L j \f~K and absorb 
j/\/K as a, while the overall energy scale J is normal- 
ized to one. Because of the large value of 7, a small 
displacement can lead to a substantial difference in the 
exchange energy. Only the adiabatic limit of the phonon 
dynamics is considered as we are ignoring the kinetic en- 
ergy. The model is defined for the two-dimensional tri- 
angular lattice with the nearest-neighbor exchange inter- 
action (ij) only. A similar model on a linear chain and 
on a square lattice has been studied earlier [12]. We look 
for the ground state of the above Hamiltonian within the 
Schwinger boson mean-field theory (SBMFT) [13] . 

One re-writes the above Hamiltonian using a pair of 
Schwinger bosons (bn,bi 2 ) defined at each site, 

Si ■ Sj — > —-AfjAij, (2) 

where A^ = bnbj 2 — h 2 bji, and decouple it by introduc- 
ing the order parameter Qij = (Aij) = —Qji. Non-zero 
\Qij\ implies finite amplitude of a singlet bond for the 
(ij) link. The resulting mean-field Hamiltonian is 

h mf = -\ E E a*ai^2 - \ E E 

+ ]T ^(6+6*1 + b+b i2 25) + \ ]T ul (3) 

i i 

where we have introduced = JijQij, and j e i to in- 
dicate all nearest neighbors of i. Lagrange multipliers Aj 
are fixed by the constraint that 8(Hmf) /dXi = which 
also amounts to requiring (b^pn + b~l 2 bi 2 ) = 25 at each 
site. Energy minimization requires that Ui obeys[10, 12] 

Ui = a^e ij (S i ■ Sj) = ~^eij\Qij\ 2 ■ (4) 

jei jei 

We proceed to solve the mean-field Hamiltonian as- 
suming (i) homogeneous pairing and (ii) inhomogeneous 



pairing of the type described in Fig. 1, and compare their 
energies. 

For the homogeneous case, the singlet pair amplitudes 
are taken as Qi,i +ei = Qi,i+e 2 = Qi,i+e 3 = iQ and A 4 = 
A. The three unit vectors are defined by ei = (1,0), 
e 2 = (-1/2.V5/2), e 3 = (-1/2,-^3/2). Passing to 
momentum space, the Hamiltonian becomes 




where = QY^ a =i sm (^ 1 e Q ). The equation of motion 
for the boson operators can be derived from Eq. (5) and 
utilizing the boson commutation algebra, 




which has a pair of eigenvalues ±Ek, Ek = \/X 2 — A 2 .. 
Boson occupation numbers rii = (bf ± bn + bf 2 bi 2 ), and the 
pairing amplitudes are worked out 

k a=l 

The fc-sum is restricted to the first Brillouin zone of the 
triangular lattice, and Nk are the number of distinct k- 
points. Upon solving Eq. (7) we obtain Q = \/3S + 
0.1161 and the ground-state energy per site 

Els. --^Q 2 --^(V3S + 0.1161) 2 . (8) 

There is no energy gain from magneto-elastic coupling 
since lattice displacement is zero for the uniform bond 
amplitudes. 

Next, we consider the situation where two distinct 
bond amplitudes, Qi and Q 2 , corresponding to dotted 
and thick lines in Fig. 1, exist. Non-zero lattice displace- 
ment develops with magnitude |itj| = (V3a/2)(Q 2 — <3i)> 
and the directions consistent with the uniform expan- 
sion of those triangles in Fig. 1 lacking a valence bond. 
Exchange energy is similarly modulated, J\ = 1 — 
(3a 2 /2)(Q 2 - Q 2 ), J 2 = 1 + (3a 2 /4)(Q 2 _ Qfy and fi _ 
nally, Ai = J 1 Q 1 and A 2 = J 2 Q 2 . 

Corresponding equation of motion of the operators 
can be worked out. Six-dimensional eigenvectors Vfc = 
(u A k u B k u C k v A k v B k v C k) T satisfy an equation of mo- 
tion Ekipk = Mkipk where 
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(9) 



2?7ifc 


= A ie jfc 


ei 


+ A 2 e lk 


e2 


+ A 2 e 4fc 


e3 




2??2fe 


= A 2 e jfc 
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= A 2 e tk 


ei 
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(10) 



Each eigenvector is normalized to obey (wAkl 2 + |wsfc| 2 + 

\u C k? ~ \vAk? ~ \v B k\ 2 ~ \v C k\ 2 = 1- 

While the problem is not tractable analytically, it can 
be shown that the (extreme) limit of Aj. = yields an 
analytical answer. We call it the "kagome limit" , where 
the non-zero A's form a network topologically equivalent 
to the kagome lattice. Using the linear Taylor expansion 
of exchange energy J\ and J 2 , the ground state energy is 
given by 



E 



1 



g.s. 



Qi 



Qi 



-a 2 {Ql-Q 



2x2 
l) i 



(11) 



for the lattice-distorted situation. In the kagome limit 
the self-consistency equations give Q\ — V3S — 0.2128, 
andQ 2 = v / 35+0.1823. Finally, in order to make J\ = 0, 
a must be equal to 0.4946. In comparing Eq. (8) with 
Eq. (11) with 5 = 2, the state with modulated bond is 
found to have a lower energy. 

The stability of the bond-modulated state can be ar- 
gued on a more general basis. Restoring the various en- 
ergy scales, the ground state energies of homogeneous 
and modulated bonds are given by 



E, 



--^JoQ 2 ! 



E™ = 



-\jiQl-J2Ql + ^u 2 . (12) 



As our analysis in the kagome limit showed, Q depends 
slowly on the amount of distortion u, and all three Q 
values are given by Q w Qi « Q 2 « \/35. Much more 
dramatic dependence on u arises from J\ and J 2 , which 
reads J\ = Jo/(l + \/3u) 7 and J 2 = Jo/(l — V3w/2) 7 as- 
suming small u. Using Q — Q\ = Q 2 = \/3S = 2 v / 3, the 
energy difference Eg s — E g n s is shown to be a monotone 
decreasing function of u for Kj Jq < 277(7 + l)i indicat- 
ing a lattice instability. The quartic term in the potential 
energy for u not included in Eq. (1) will restore the sta- 
bility at a finite u and favor the lattice reconstruction 
indicated in Fig. 1. 

The mean-field ground state features Bose conden- 
sation for both uniform and non-uniform states which 
are the hallmarks of magnetic order. For k — k* at 



FIG. 2: Static structure factor S(k) = (S-k ■ Sk) for the 
|AKLT) state in Eq. (13). Solid line is S(k, 0), and the dashed 
line is 5(0, k). MC calculations were performed on N = 90 x 
90 punctured triangular lattice. 



which the Bose condensation occurs, the ground state 
is described by a form |LRO) = (7^)^*10) {N k , = 
number of bosons in the condensed state). For other 
k values the ground state is defined by 7fe|MF) = 0, 
where 7^ are the quasiparticle operators diagonalizing 
the mean- field Hamiltonian, Eq. (3). The uncon- 
densed part, written in real space, is of the general 
form, |MF) = exp <jj - 6+6+)} |0), rep- 

resenting a gapped spin liquid[13, 15]. The sensitiv- 
ity of the ground state on a particular spin quantum 
number 5* arises when we make a projection to a fixed 
S: |VBS) = P<?|MF). The full ground state is a di- 
rect product of the condensed and the uncondensed parts, 
|GS) = |LRO) <8 |VBS). 

With 5 = 2, the projected ground state consists of 
a large number of valence bond configurations all shar- 
ing the property that each site has four valence bonds 
emanating form it. Among those states, the VBS state 
depicted in Fig. 1 is the one most consistent with the 
kagome limit. Such a state, in turn, is described by the 
AKLT state: 



|AKLT) 



'(2 



tfi)|0>. 



(13) 



(ij)' runs over the four valence bonds depicted in Fig. 1. 

Although the equivalence of |AKLT) with |VBS) S=2 
cannot be established rigorously, it is still reasonable 
to expect that the physical properties embodied in the 
I AKLT) state will reflect those of the |VBS) S=2 [16]. 

The quasiparticle spectrum obtained from diagonaliz- 
ing Eq. (9) reduces to the familiar spin wave modes[14]. 
The excitation energies for the AKLT ground state given 
above can be worked out in a manner parallel to the cal- 
culation for the gapped spin liquid in one-dimensional, 
5=1 chain[6]. We adopt the method of single mode 
approximation (SMA), and write the excitation energy 
A(fc): 



A(k) 



([5_ fc -,[g,S fc ]) 
2(5_ fc • S k ) 

1 - (1/3) ELi COf # ' e c 



rj 1 k ' ( Tj 



(14) 
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The averages (•••) are to be taken with respect to 
|AKLT) state given in Eq. (13). The nearest-neighbor 
spin-spin correlation —B\, as well as all other averages 
(Si ■ Sj) can be evaluated using the standard Monte Carlo 
(MC) method [15, 17]. The effective coordination number 
on the kagome lattice is z' = 4; the number of sites of 
the triangular lattice is N. 

The structure factor <S(fc) = (S-k • Sk) vanishes in the 
k = limit due to the singlet character of the AKLT 
state. The next term in the small- fc expansion gives 
S{k) « Tfc 2 , where -T = X^.fo - r,) 2 (^ • Sj). 

The numerator in Eq. (14) also vanishes as k 2 , giving 
the k = gap of magnitude A(k = 0) = 8 J x (Bi/T). 

The AKLT state has only short-range order as mani- 
fested by an exponential decay of the spin-spin correla- 
tion, \(Si ■ Sj)\ ~ e" 1 ' 467 !^'!, in our MC simulation. In 
Fig. 2 we present results for the static structure factor. 
The k = gap obtained from S(k) is A(0) ~ 0.5 x (6J). 

Figure 3 displays excitation spectrum A(fc) for |VBS) 
according to Eq. (14), and the standard spin wave in the 
uniform |LRO). The new spectrum A(fc) is characterized 
by a broad minimum (dark region) that exist around the 
zone boundary of the Brillouin zone, and in principle 
should be observable in the neutron scattering. The spin 
waves and the new excitation mode we find reflect the 
two kinds of broke symmetries found in the ground state. 
Interaction between the two modes should exist, which 
requires going beyond the quadratic Hamiltonian (3) to 
include fluctuation terms. 

In summary, we propose that a VBS state can be sta- 
bilized in two-dimensional antifcrromagnetic spins with 
a mismatch of 2S and z, when a sufficiently strong cou- 
pling to the lattice deformation exists. As an example, 
we analyze S = 2, z = 6 (triangular) model in the 
Schwinger boson mean-field theory and show that the 
energy is lower for the lattice-deformed state. Corre- 
sponding gapped spin excitation energies are calculated 
within the single-mode approximation. 
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